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The problem of geometric phase for an open quantum system is reinvestigated in a unifying approach. Two of 
existing methods to define geometric phase, one by Uhlmann's approach and the other by kinematic approach, 
which have been considered to be distinct, are shown to be related in this framework. The method is based 
upon purification of a density matrix by its uniform decomposition and a generalization of the parallel transport 
condition obtained from this decomposition. It is shown that the generalized parallel transport condition can 
be satisfied when Uhlmann's condition holds. However, it does not mean that all solutions of the generalized 
parallel transport condition are compatible with those of Uhlmann's one. It is also shown how to recover the 
earlier known definitions of geometric phase as well as how to generalize them when degeneracy exists and 
varies in time. 
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The concept of geometric phase was originally introduced 
by Pancharatnam in the classical context of comparing two 
polarized light beams through their interference |1]. Later, 
Berry pointed out to its importance even in quantum systems 
undergoing a cyclic adiabatic evolution |2]. After that, this 
important notion has been subject of interest in many differ- 
ent aspects, which has led to many different generalizations 
and applications 01 01- Of course in general cases to retain 
purely geometrical nature of the phase one has to put some 
constrains, namely parallel transport (PT) conditions. In this 
manner, geometric phase is a feature which only depends on 
geometry of the path traversed by the system in its motion 
during evolution. 

It is also worth noting that an important source of the re- 
newed interest in geometric phases is their relevance to geo- 
metric quantum computation and holonomic quantum compu- 
tation 1 5]. Indeed, it is known that quantum logic gates can be 
implemented only by using the concept of geometric phases. 
It is believed that the purely geometric nature of this phase 
makes such computations intrinsically fault tolerant and ro- 
bust against noise |6|. 

A pure state is merely an idealization and in real experi- 
ments a description of the system in terms of mixed states is 
usually required. This point accounts for attempts towards ex- 
tending the concept of geometric phase to mixed states. In 
fact, Uhlmann was the first in tackling the problem through 
the mathematical approach of purification of mixed states |7]. 
This method is rather general in that it is independent of the 
type of evolution of the system. Next, Sjoqvist et al. put 
forward a quantum interferometric based definition for geo- 
metric phase of nondegenerate density matrices undergoing a 
unitary evolution 1 8]. Later, Singh et al. proposed a kinematic 
description and extended the results to the case of degenerate 
mixed states [9]. It must be mentioned that there also exists 
another, differential geometric, approach to define geometric 
phase for mixed states undergoing a unitary evolution fioll . 
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In this approach, mixed state geometric phase appears as an 
immediate and direct generalization of the pure state case. 

Indeed, in the case of environmental effects such as deco- 
herence, one has to consider nonunitary evolutions of mixed 
states. Some generalizations in this direction have been ad- 
dressed in Refs. . The 
proposition in Ref. 1 12] for completely positive maps in spite 
of being operationally well-defined depends on the specific 
Kraus representation for the map. In Refs. fTMfl^l . the prob- 
lem of geometric phase of an initially pure open quantum sys- 
tem, based on the standard definition of pure state geometric 
phase, has been addressed through the quantum jump method. 
A more recent effort is based on a kinematic approach, with 
no a priori assumption about dynamics of the system 0. 
However, most of these different definitions do not agree with 
each other. In fact, Uhlmann's method even in the case of 
unitary evolution does not agree with the interferometric defi- 
nition l2lll22Tl . The source of such disagreement is known to 
be the use of different types of PT conditions. Hence, it has 
been argued I2ll I22I l23ll that these approaches are not gen- 
erally equivalent and one cannot obtain one from the other. 
Therefore, it could be desirable to find a more unified ap- 
proach which can bring together the previous general ideas. 
Recently, in the unitary evolution case it has been argued that 
using (nonorthogonal) decompositions different from spectral 
decomposition can make it possible to unify the kinematic and 
Uhlmann's approaches 12411 . In this framework, a suitable no- 
tion of PT condition of the mixed state is based on the PT 
condition of the vectors constituting this decomposition. 

In this paper, we shall use a rather similar mechanism plus 
uniform decomposition of density matrices, and propose a 
generalized kinematic approach for geometric phase of mixed 
states under an arbitrary nonunitary evolution. This approach 
vividly shows how it is possible to merge Uhlmann's approach 
and kinematic approach. It is also shown how to recover the 
earlier definitions of geometric phase from this more general 
approach. In addition, it is shown that the approach can be 
easily modified to include the more general case of degener- 
ate mixed states. This investigation may as well be useful in 
the study of robustness of geometric phases against decoher- 
ence 1 2511. 
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Let us suppose that the density matrix of our sys- 
tem of interest (with the Hilbert space TL S ) is g(t) — 
Y,k=iPk( t )\ w k{t))(w k (t)\, in which Pk(t)'s (\w k (t))'s) are 
considered to be its eigenvalues (normalized eigenvectors). In 
a general evolution both pk and \wk) are subject to change 
in time. For simplicity of our discussion, we in the se- 
quel assume that the rank of this matrix is constant at all 
instants, and even more the matrix is nondegenerate. In 
the case of unitary evolution, we have Pk{t) = Pfc(O) and 
\ w k(t)} — U(t)\wk{0)), where U(t) is the unitary evolution 
operator. However, when evolution is nonunitary, the eigen- 
values, p/,, can also vary in time. Thus, generally U(t) = 
J2k \ w k(t))(wk(0) \ does not encompass the whole dynamical 
information. In fact, in such cases, to obtain g(t) one often has 
to resort to some approximative methods in the theory of open 
quantum systems 1 19], such as the Lindblad equation 1 20] . 

Since in our construction we use Uhlmann's PT condi- 
tion |7|] we need to recall it briefly. Uhlmann's approach 
is based upon the standard purification w(t), where g(t) = 
w(t)w^(t), for density matrices. In other words, w can be 
considered as a purification of g in the larger Hilbert space of 
Hilbert-Schmidt operators with scalar product (w(t), w(t')) = 
ir(w^ (t)w(t')) such that = g. It is clear that w(t) — 
\J g{t)V(t) is an acceptable purification of g for any unitary 
V(t). For a special purification where each \(w(t),w(t'))\ 
is constrained to its maximum value Uhlmann has defined 
the geometric phase associated to the evolution from g(0) to 
g(r) as 7 fl (r) = arg((w(0), w(r))), where the PT condition 
w^(t)w(t) — w^(t)w(t) has to be satisfied. 

Let us also briefly review the construction of the geometric 
phase in Ref. I17ll . Consider a purification for the density 
matrix g(t) as 

\y{t)) S a=EkVPkTt)\Mt)}s®Wk)a, *G[0,T]. (1) 

Now after imposing the PT condition; (wk(t)\^\wk(t)) = 0, 
the geometric phase, defined a la Pancharatnam 
d; 7(1") = arg((# (0)|W(-r))), reads as 7ff (r) = 
arg(Ef=i VPfc(0)PfcW(^(0)|«; fc (r))e-/o T <«'*W|ti-*(t)>dt). 
Indeed, by using the PT condition one fixes the general form 
of the unitary operators which like U(t) can run system's 
dynamics. As is clear, in this method purification of mixed 
state of the system is done based on its spectral decomposition 
and the PT condition is considered to be the PT condition 
of all the vectors constituting this (spectral) decomposition. 
We know that a purification as in Eq. 0, is only one of the 
possible purifications that can give rise to the correct mixed 
state of the system. So, one has the freedom to choose other 
decompositions and study the problem of geometric phase 
with respect to them. In the sequel, we follow such a strategy 
and look for a specific purification in which all normalized 
terms can be treated in a naturally uniform manner, unlike 
Eq. where the contribution of the fc-th normalized term is 
the time dependent variable \/pk(t)- In other words, instead 
of starting from the spectral decomposition of a density 
matrix which is the usual starting point of purification based 
approaches, we start with another decomposition which can 
result to the mentioned uniformity. In order to do so, we need 



the next two important theorems on different decompositions 
of a density mat rix g . 

Theorem 1 1 26] : Let g has the spectral ensemble 
{Pk, \u>k) }■ Then {qi, \xi)} is another ensemble for it iff there 
exists a unitary matrix U = (JAu ) such that 

= J2kVPk^k\wk)- (2) 

Theorem 2 ll27ll : Let {qi} be a probability distribution. 
Then there exist normalized quantum states such that 

Q = J2i qi\xi)(xi\, iff qis majorized by p. 
An immediate corollary of Theorem 2 is the existence of a 
uniform ensemble for any density matrix. Therefore, there 
exist normalized pure states \x%), . . . , \xj\f) such that g is an 
equal mixture of these states with probability 1 /Af (Af > N), 

i.e. g = jf Ym=i \ x l){ x l\- F° r tr| e rest of discussion we as- 
sume that Af = N. Now, let us see how this uniform de- 
composition is related to the spectral decomposition. By us- 
ing Theorem 1, we have -^\xk) = Ya=i ^/piU M \wi). It 
is easy to see that if one chooses an N x N Fourier ma- 
trix (corresponding to discrete Fourier transformations 12811 ') 
U k i = 7^ e_27 ™^ (M = 0, . . . , N - 1), and momentar- 
ily run all indices from to N — 1 (rather than 1 to N) this 
equation is satisfied. Then, by using a Fourier matrix one can 
find a uniform ensemble for any density matrix. If we define 
C(t) = J^k yPk(J)\wk(0)}(wk{0)\ and use the definition of 
U(t), we can rewrite \xk(t)) in the following matrix form 

\x k (t))=VNU(t)C(t)U\w k (0)). (3) 

Now we show that the above mentioned uniform decompo- 
sition is useful in our discussion of geometric phase. Consider 
the following pure state of the combined system sa 

mt)) sa = J= J2 M)) s ® V(t)\a k ) a , (4) 

where V(t) is the unitary evolution of the |ofc)'s. This state 
is a legitimate purification of the density matrix g(t) of the 
system; g(t) = tr a (|$(t)) sa ($(*)|). If V{t) = I, since 
(xk\xk) — 1 and all \xk(t)} vectors enter with equal and 
constant probability of -k in the decomposition of the den- 
sity matrix, it seems natural to consider our (generalized) PT 
conditions in the form of 

(ar ik (t)|^|a: fc (t))=0 ) k = l,...,N, (5) 

that is, a density matrix undergoes a PT condition when all 
of the vectors in its uniform decomposition do so. Here a 
point is in order. It must be mentioned that, except the pure 
state case, this PT condition is generally different from the 
one considered in earlier literature |8|,[l7|] 1 . 



Tong et al.'s PT condition llTll imposes a constraint on the 

form of parallel transported evolution operator, U^(t), as 

(iu fc (0)|t7llti/ll|«; fc (0)) = 0, k = 1, ...,7V, whereas Eq. HJ 
results into {w k (0)\W CU^U^CU\w k (0)) =0, k = l,...,N. 
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In general, in the purification (|4} ancillary vectors could 
also vary in time, and we have to find a natural picture for geo- 
metric phase in this case. Let us first remind a simple and use- 
ful p roperty of Schmidt decomposition of bipartite pure states 
M- If |*}„& = Ek c Ma\b k )b, then (U®V)\$) ab = 
(UCV T <8 I)E fe l dk)a\bk)b, where C is a diagonal matrix 
in the {|afc)} basis defined as C = i~2 k c k \a k )(a k \ and V = 
Efcfc' (t>k\V\bk>) \a k ) («fc' I- Here, for notational purposes, we 
omit T sign of V T . Now, noting this property and assum- 
ing that the basis vectors of the ancillary Hilbert space are 
{|wfc(0))}, one can rewrite Eq. (0} as 

!*(*)>-. = Ek !«*(<))» ®l«*(o))«» ( 6 ) 

where \x k (t)) = U{t)C(t)UV(t)\w k (0)). This purifi- 
cation now results into the nonorthogonal decomposition 
g(t) = J2k \%k(t)) (xk(t)\ for the density matrix. Unlike the 
{|£fc(i))} decomposition, now for a general V, (x k (t)\x k (t)) 
is not time independent and, as well, is not equal for all fc's. 
However, if we consider the normalized vectors \x(t)) — 
J - fc ^n it still looks natural to consider our generalized PT 
condition to be in the following form 



(xk(t)\-\x k (t)) 



0. 



(7) 



In terms of \x k (t)) vectors this is equal to {xk(t)\^\x k (t)) = 
T}-jt((xk(t)\xk(t)}), or equivalently in more detail it is 

{w k {Q)\V^U^CU^UCUV + V^U^CCUV + V^C 2 U 

Id., ..... (8) 



xVK(O)) 



2dt 



((w k (0)\VWC 2 UV\w k (0))). 



Now let us see what is the form of Uhlmann's PT condition. 
We note that w(t) operator reads as w(t) = U(t)C(t)UV{t). 
Hence, the explicit form of Uhlmann's PT condition is 

vWctfucuv + vWccuv + V ] U^C 2 UV 

(9) 

= vWctfucuv + vWccuv + vWc 2 uv. 

As is seen LHS of this equation is exactly the expression 
within bra-ket of the PT condition (jHJl. If sandwiched between 
(wfc(0)| and \w k (Q)), Eq. (|9) gives rise to 



LHS of JgJ = i(w fe (0)|LHS + RHS of @J|l£>fc(0)) 



1 d 

2dt 



((w k (0)\vWC 2 UV\w k (0))) 



(10) 



This is what we wanted to show; by using Uhlmann's PT 
condition the generalized PT conditions (|8jl are also satisfied. 
However, it must be noted that generally number of equations 
of the two PT conditions are not equal. In other words, Eq. (|9} 
is a matrix equation which constitutes N 2 different equations 
(for V) though Eq. is just a set of N equations. This sim- 
ply means that there might be solutions of Eq. (|8j that are not 
solutions of Eq. (J9ji- If it is assumed that V(i) = e^ lH ^\ then 
solution of Eq. l|9} is as follows |7] 



-iH(t) 



2£V| uv (0))(u; fc (0)|M 

kk' 

x / dt {w k '{t )\w k (t ))- 



(11) 



Now it is easy to show that Eq. (|8j can have solutions other 
than (II It . For example, if we suppose that [UV, C] = 0, and 
UV = E fe e" 8l(fl K(0)) (w k (0)\, then Eq. © gives 



h(t) = -i jo'df '{w k (t')\w k (t')}, 



(12) 



which does not generate a V{t) compatible with (II li . This 
comes from the fact that to satisfy Eq. (|SJl we only need to 
have the diagonal terms of Uhlmann's PT condition, whereas 
off-diagonal terms of this equation may put extra constraints 
that are redundant for validity of Eq. (|8j- 

Now geometric phase can be simply defined a la Pancharat- 
nam as 

lg {t) =arg(($(0)|*(t)» =arg(£ fe ^)e l7fc(t) ), (13) 

where (x k (0)\x k (t)) = i/ fe (t)e^( t ), i.e. v k {t) ( 7fc (t)) is the 
visibility (geometric phase) of the fc-th component of 
The explicit form is obtained by insertion of the definition of 
\x k (t)) in this equation, which gives 



7 fl (t) =arg(J] ^p k >{Q)p k {t){w k .{Q)\w k {t)) 



kk' 



(14) 



x (w k (0)\UV(t)W\w k ,(Q))). 



This equation shows that geometric phase, as described here 
to be combined with Uhlmann's definition, generally retains a 
memory of the evolution of both system and the ancilla, that 
is, it is a general property of the whole system which depends 
on the history of the system as well as the history of the ancilla 
entangled with it 12 ill . 

In the remainder, we investigate how the earlier defini- 
tions of geometric phase 0, fl7ll can be obtained from the 
present framework as special cases. If we confine ourselves 
to a restriction of the solution of {^l for V(i), such that 
Vfcfe' (t) = Vkk' (t)S k k' and has the property 



[V,WC 2 U] = 0, 



(15) 



(i 



p k '{t')+p k {t')' 



or equivalently V(t) = J2k£_ ilk(t) ^\ w k(0))(w k (0)\U, 
where l(t) is defined as in Eq. dl2> . then explicit form of j g 
becomes 

lg(t) - arg(£ fc ^/ Pk (0)pk (t)(w k (0)\w k (*))e-«*«), (16) 

as in Ref. fnll . Thus, in the context of the discussion of 
Ref. 0, it can be said that the physical role of the commu- 
tation relation dl5> appears like removing memory effects of 
ancilla's evolution from geometric phase. 

Let us end by mentioning some remarks on the initial 
assumptions of the approach while our stress is still on 
derivation of earlier results and their possible generaliza- 
tions. Based upon Theorem 2, it is seen that one can al- 
ways choose Af, number of the vectors in uniform decom- 
position, such that Af > N. For example, we can assume 
that Af = dim(Tt s )- Now we show how the whole frame- 
work can be modified in the degenerate case. Consider the 
evolution for the density matrix of the system from g(0) 

to g(f) = Ef=iEZliPfc(*)K(*)>K(*)l. where Pk (t), 
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k = 1, . . . , N, are the n^-fold degenerate eigenvalues of g(t), 
and \w^(t)), fi — 1, . . . , nf., are considered the correspond- 
ing eigenvectors. In this case, the pure state of the total sys- 
tern is | *(<)>„, = £f=i YTjLi !*£(*)>, ® K(0))„, where 
|$t(t)) is defined as in Eq. lUJi in which \wk(0)) is replaced 
by \w%(0)). Then, one notes that 

($(o)|$(t))= £ vSM)W(o)l<W) 

kk'nn' (17) 

x (^(0)|WV(i)WtK'(l))). 

which is determined when all elements of V(f) are known. 
Now we choose our PT condition in this general case as 

(5fc(*)l^$f(*)>=0, /i.M^l,-.,"*. (18) 

It can be checked that this PT condition can also be satisfied 
by assuming Uhlmann's PT condition, Eq. (|9}. In this case, it 
is easily seen that the most general form for V which satisfies 
Eq. O is as follows 

v(t) = E w oT'W z^K(o))«'(o)|w. (i9) 

After some algebra and using the commuta- 
tion relation H5\ . it is obtained that (t) = 
(w^(0)\Pe-^ u Ht')u(t'yt'\ w ^' (p)) t where p denotes 
path ordering. After inserting this relation back into Eq. d!7i . 
non-Abelian factors show up in the geometric phase. 

In general, when degeneracies vary in time, a level- 
crossing like behavior can occur. In this situation, in the dis- 
cussion of differentiability of the eigenvalues (and eigenvec- 



tors) the notion of ordering of the eigenvalues becomes impor- 
tant. For example, it can happen that the natural ordering as 
pi(t) > ■ ■ ■ > pN(t) (for all t) destroys differentiability, thus, 
one has to seek for some ordering which respects it 12911 . If 
such an ordering can be found, then the operator U(t), eigen- 
values, and eigenvectors are still well-defined differentiable 
functions and our approach may be generalized as well. 

In summary, the notion of geometric phase of a mixed state 
undergoing nonunitary evolution has been investigated in a 
unifying picture in which two of the previous general defi- 
nitions, Uhlmann's definition and kinematic approach, have 
been related to each other. In this formalism, we have used 
the idea of purification of state of a system by uniform de- 
composition of its density matrix rather than the spectral one, 
and by attaching a time varying ancilla to it. Then, as a natu- 
ral choice for parallel transport condition, we have considered 
that a mixed state is undergoing parallel transport condition 
when all the (normalized) vectors of its corresponding purifi- 
cation are subject to this condition. This generalized parallel 
transport condition is different from the ones defined previ- 
ously in the literature. It has been shown that the new condi- 
tions are satisfied when Uhlmann's condition holds. However, 
because of different numbers of equations in the two parallel 
transport conditions, the generalized parallel conditions are 
only diagonal equations of Uhlmann's condition. Finally, it 
has been shown how to recover earlier definitions of geomet- 
ric phase of a mixed state. Extension of the method to the 
more general cases of degenerate density matrices with time 
varying degeneracies have also been discussed. 

We thank N. Paunkovic for useful discussions. This work 
was supported by EU project TOPQIP under Contract No. 
IST-2001-39215. 
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